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INVESTIGATION OF TEMPERATURE MEASUREMENTS 
IN 300° to 1100° K LOW- DENSITY AIR USING 
AN ELECTRON BEAM PROBE* 

By William W. Hunter, Jr. 

Langley Research Center 

SUMMARY 

Laboratory measurements of rotational and vibrational temperatures have been 
performed over a range from 300° to 1100° K in static low-density air with an electron 
beam probe. The tests were conducted in a unique test chamber which permitted the 
selection of various combinations of test pressures and temperatures. 

Rotational-temperature measurements were made from the (0,0) band of the first 
negative system of the nitrogen constituent of air. Some tests were also performed by 
using the (0,1) band of the first negative system of nitrogen. No substantial differences 
were noted between the results of the (0,1) band tests and those of the (0,0) band tests. 

The accuracy of the measured rotational temperatures was found to vary with gas 
temperature and gas number density. A difference existed between the measured rota- 
tional temperature and the gas temperature. This difference increased with increasing 
gas temperature. With the inclusion of an experimentally determined correction for the 
gas -number -density effect, the rotational-temperature measurements were made with an 
accuracy that varied from approximately 0 to -6 percent with increasing gas temperature. 

Vibrational-temperature measurements were made by using the intensity ratios of 
the (0,1) and (1,2) bands of the first negative system of nitrogen. The measurements 
were found to be within ±20 percent of the gas temperature. 

INTRODUCTION 

The present investigation was undertaken as part of a continuing instrument develop- 
ment program of the Langley Research Center. This particular effort was directed 

Part of the material presented in this report is included in a thesis entitled 
"Measurement of Rotational and Vibrational Temperatures Over a 300° K to 

1100° K Range Using a High-Energy Electron Beam" submitted in partial fulfillment of the 
requirements for the degree of Master of Arts in Physics, The College of William and 
Mary in Virginia, August 1965. 


towards the verification of a technique for measuring free- stream temperatures of low- 
density hypersonic wind tunnels which use air or nitrogen. 

The technique consists of passing a beam of high-energy (10- to 30-kV) electrons 
through a low-density gas. The electrons have inelastic collisions with gas molecules 
and these excited ionized molecules provide fluorescence. A spectral analysis of the 
fluorescence provides a means for determining molecular gas rotational and vibrational 
temperatures. 

The initial investigator, E. P. Muntz (ref. 1), experimentally verified his theoreti- 
cal model for the electron beam probe technique for two temperatures, =300° K and 
373° K. These experiments were conducted in a low-velocity (0.5-m/sec) flow of nitro- 
gen gas, and the data were obtained from the (0,0) band of the first negative system of 
nitrogen. Since the original work, others (refs. 2, 3, 4, and 5) have used the technique 
for diagnostic wind-tunnel measurements. In reference 4, the calculated tunnel param- 
eters are compared with the measured rotational temperatures as a test of the theory for 
temperatures below 300° K. 

The purpose of the present work was to verify experimentally the theory for 
rotational- and vibrational -temperature measurements for a range of temperatures from 
300° to 1100° K under controlled conditions in the laboratory. These tests were conducted 
in a static test gas, air. The rotational-temperature measurements were made from data 
obtained from the (0,0) band of nitrogen. Rotational-temperature measurements were 
also made from data obtained from the (0,1) band of the first negative system of nitrogen 
for comparison with results obtained from the (0,0) band. The vibrational-temperature 
measurements were made from data obtained from the (0,1) and (1,2) bands. Tests were 
conducted in nitrogen as well as in air to determine whether there were any significant 
differences in the resulting temperature measurements. 

SYMBOLS 

(A) term defined by equation (25) 

Ajt’k'^ transition probability for emission between rotational energy states K 
and K *2 (see eq. (23)) 

A nm transition probability for emission between states n and m 

A v Vo transition probability for emission between vibrational energy states v' 
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rotational constant related to vibrational level v'^ = 0 

+ O + 

electronic wave function for N2B^S U 

excitation function which described electron-molecular excitation process 
ratio of excitation function to depopulation rate of N^B^S* state 
speed of light 

characteristic energy of rotational energy level 
characteristic energy of vibrational energy level v T ^ 
electron charge 

initial wave function of primary electron 
final wave function of primary electron 
wave function of secondary electron 
rotational term (ref. 9) 
defined by equation (30) 
vibrational term (ref. 9) 

Planck T s constant 

intensity of emission for transitions between K f and rotational energy 

levels 

intensity of emission for transitions between n and m states 
reference intensity of emission 

intensity of emission for transitions between v T and v^ vibrational energy 
levels 
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quantum number of total angular momentum 
rotational quantum number 
change of rotational quantum number 
Boltzmann's constant 

quantum number of a component of total angular momentum 

steady- state number density population of a rotational energy level 

number density population of state n 

steady-state number density population of N 2 xlXg 

steady-state number density population of a vibrational energy level 

neutral nitrogen 

ionized nitrogen 

excited ion state of N 2 

ground state of N 2 

ground state of N 2 

gas number density 

reference gas number density, 3 x 10” ^ cm -2 
relative rotational line strength for excitation 

relative rotational line strength for excitation for R branch and P branch, 
respectively 

relative rotational line strength for emission 
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P y i v < i band strength 

Q r rotational partition function 

Q v vibrational partition function 

q v > v 't Franck-Condon factor 

q momentum transfer term 

R depopulation rate of N*jB S u 

R^j electronic transition moment for electronic states i and j 

r n nuclei separation 

r^ distance between primary and ith orbital electrons 

r^2 distance between primary and secondary electrons 

r^ position vector of ith orbital electron 

r j position vector of primary electron with respect to point of observation 

?2 position vector of secondary electron with respect to point of observation 

r 2 position vector of secondary electron with respect to molecular axis 

X£ T 

Sjr" Honl-London factor (ref. 9) 

T gas temperature 

T r rotational temperature 

T v vibrational temperature 

v vibrational energy state 

W constant defined in equation (26) 
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X 


constant in equation (23) 


X^Sg electronic wave function for N 2 X^Eg 

Y constant defined in equation (29) 

Z constant in equations (1) and (31) 

9,x,<P Euler angles 

A quantum number of resultant electronic orbital angular momentum 

v wave number 

reference wave number 
wave number of (n,m) transition 
wave number of (v^v'rjj transition 


^o 


'nm 




°u 2s 




angular characteristics of an orbital electron of N 2 XlEg (ref. 9) 
electronic- state wave function 


i^evjAM molecular wave function 

p jam rotational- state wave function 

i p v vibrational-state wave function 

Superscripts: 

' excited ionized state of nitrogen, 


ground state of neutral nitrogen or ionized nitrogen 


Subscripts: 

a vibrational energy state a of excited ionized nitrogen 

b vibrational energy state b of excited ionized nitrogen 
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1 ground state of neutral nitrogen, N 2 X x £g 

2 ground state of ionized nitrogen, ^X^s* 

EXPERIMENTAL SYSTEM 
Test Chamber and Control System 

The test chamber and control system, shown in figure 1, was designed to provide 
flexibility in temperature and vacuum test conditions. Temperature and pressure oper 
ating ranges are approximately 300° to 1100° K and 6.7 to 133.3 x 10“^ N/m2. 



Figure 1.- Experimental apparatus. L-66-7099.1 
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The major component of the system is the test chamber (fig. 2), which consists of 
three concentric cylinders. The outer cylinder is a stainless-steel water-cooled jacket 
and is fitted with vacuum-tight water-cooled top and bottom covers. Each cover is fitted 
with large flanges, attached to extensions, for mounting test hardware. Three 7.6-cm- 
diameter optical-grade-quartz windows are located in the outer cylinder wall. The next 
concentric cylinder consists of a helically wound nickel-ribbon heating element and is 
attached to ceramic supporting rods. Electrical connections are made to copper 
electrodes which extend through the outer cylinder. The inner cylinder is a 20-cm- 
diameter, 38-cm-long stainless -steel electrostatic shield and is grounded together with 



Figure 2.- Test chamber and electron gun. 



the outer cylinder to prevent charge buildup on the walls. Also, the inner cylinder pro- 
vides a more uniform heating surface for the test gas than would be provided by the rib- 
bon heating element. The inner cylinder is equipped with end covers which have openings 
for passage of the electron beam. A 3.8- by 1.3-cm slit opening is provided in the cylin- 
der wall and is located in line with one of the viewing windows of the outer cylinder. 

Rectified heater current is supplied from a 440-V ac three-phase system. Temper- 
ature control is provided by coarse and fine rheostats. Temperature is regulated within 
±1 percent of the preset value by an on-off automatic pyrometer. An operating tempera- 
ture of 1100° K is obtained for a heating- element voltage of 35 V dc and a current of 
50 amperes. 

A 35-psi (241-kN/m2) water-cooling system provides a heat sink for the outer 
chamber wall and covers and cools the diffusion pump. An interlock system prevents 
operation of the heating system and diffusion pump unless proper cooling flow is 
established. 

A 5-ft^/min (0.002-m3/sec) mechanical pump, 750-liter/sec (0.75-m3/sec) diffu- 
sion pump, cold trap, and necessary isolating valves comprise the vacuum pumping sys- 
tem. A variable leak valve with air dryer, in combination with the mechanical pump, was 
used to maintain the test chamber at the desired pressure. For all experiments, the 
chamber pressure was measured with a McLeod gage. 

Test-Chamber Temperature Survey 

A series of tests was conducted to establish the test- chamber gas temperature. 
First, a survey of the inner-cylinder wall was made to determine the temperature varia- 
tion of the heating surface. This survey was performed with 12 thermocouples attached 
to the inner-cylinder wall and located to give reasonable coverage. The results indicated 
that a ±20° variation existed at 1100° K. This variation decreased with decreasing 
temperature. 

The inner-cylinder ends were not in close proximity to the heating element and 
were, therefore, cooler than the side walls. The second test was made to determine the 
effects of these cooler surfaces on the test gas at the point of observation. A thermo- 
couple was located on the bottom end plate near the chamber center line. The tempera- 
ture of this point was the lowest temperature of the plate. The temperature along the 
center line was calculated from a solution of Laplace's equation in cylindrical coordinates. 
It was assumed in this solution that the temperature varied linearly along the length of 
the cylinder and was constant about the circumference. This assumption was justified 
on the basis of the data obtained from the tests described in the preceding paragraph. 
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The results of the calculations show that the effects of the inner -cylinder ends would 
lower the gas temperature at the midpoint of the cylinder no more than 4° at 1100° K. 

The final test was a direct measurement of the test gas translational temperature 
(hereafter referred to as gas temperature) at the point of observation. These measure- 
ments were made with a thermocouple whose leads were brought in along the chamber 
center line and the junction located at the point of observation. Two thermocouple sizes, 
0.1-mm and 0.5-mm diameter, were used to determine whether there were significant 
heat losses due to thermocouple conduction at the point of observation. The data show 
that conduction losses were negligible. The results are indicated in figure 3, which is a 
plot of the difference between the gas temperature, measured with the thermocouples, and 
the reference wall temperature as a function of the test-chamber inner -cylinder reference 
wall temperature. These data show that the gas temperature at the point of observation is 
lower than the reference wall temperature. This difference decreases with a decrease in 
temperature and is attributed to heat sinks provided by the observation and pumping ports. 


Gas (air) Thermocouple 

number density diometer, 

cm'3 mm 



Figure 3.- Test-chamber wall temperature as a function of difference between wall temperature and gas temperature. 

Electron Beam System 

The electron gun (fig. 4) is a conventional point cathode system which has a directly 
heated hairpin tungsten filament, independent negative grid bias for current control, and 
cathode focusing. The beam is magnetically focused and deflected upon passage into the 
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drift tube. The focused beam passes through a 1.0- mm -diameter hole in a 2.5-cm-long 
plug which is placed in the end of the drift tube. Because of the low gas conductance of 
this passage through the plug, the test chamber may be operated at pressures up to 
133 N/m.2 while the gun is maintained within acceptable pressure range, S2.7 cN/m2. 

The beam system was operated at potentials between 25 and 30 kV and currents 
between 1.0 to 1.5 mA. The beam potential and currents were held constant within ±1 per- 
cent of the selected values for all tests. 


High- 
vol tage 
cable 


Upper 

deflection 

coil 


Focus 

coil 



Lower 

deflection 

coil 



Drift 

tube 


Figure 4.- Electron gun and anode assembly. L-68-810 

Optical and Electronic Detector System 

The major component of the optical and electronic detector system is the 0.5-meter 
scanning spectrometer. This instrument has 16.0-A/mm dispersion in the first order 
and 0.2-A resolution. A 13-stage venetian-blind photomultiplier tube which has an S-13 
spectral response characteristic is mounted at the exit slit. The photomultiplier output 
is fed into an electrometer amplifier which drives a strip-chart recorder. The overall 
response of the system is approximately 1 second. Because the response of the available 
system was relatively slow, a 5-A/min scanning speed was used. 

ROTATIONAL-TEMPERATURE MEASUREMENTS 

Tests 

Rotational-temperature measurements were made in air from the (0,0) band of the 
first negative system of nitrogen at approximately 100° intervals from 300° to 1000° K. 
Two sets of tests were performed with a constant gas number density of 1 x 1 0 16 cm“3 
and the other set, with a constant gas number density of 5 X 10l5 cm _ 3. Typical spectral 
traces are shown in figure 5. For all rotational-temperature tests, the width of the spec- 
trometer entrance slits were selected so that the intensity of each line could be determined 
from the peak value of each line. 

The procedure for obtaining the rotational-temperature measurement was to deter- 
mine the relative peak value of a rotation spectral line and to enter this value for 
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Ng (0,0) band 
temp., 800 ° K 


Intensity 



into the following equation, which is derived in the appendix: 

^Ao 


Bnhc . 

-t^-K (K + 1) + Z = -2.3 log 10 
kT ■£ 


(k' + Kg + 1 




(1) 


A plot of the value of -2.3 log^g 


'k'k’sAo 


(k' + K ’2 + 1 ) \G\v/v 0 f 


- against K’(K f + 1) provides 


a datum point for each spectral line. The two curves in figures 6 and 7 are obtained by a 
least-squares fit to the data points. One curve represents the odd K' lines (strong line 
system) and the second curve represents the even K* lines (weak line system). The 


slopes of the curves are equal to 


B()hc 

~kT\ 


, and from this equality the rotational temperature 


T r is determined. The final value of the rotational temperature is determined from a 
weighted mean of the values obtained from the two curves. The weighting is in accordance 
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Figure 6.- ^(0,0) band data for 300° K experiment for T r measurements. 



Figure 7.- N^tO.O) band data for 800° K experiment for T r measurements. 


with the standard deviation of each curve with respect to its best straight-line fit to the 
data points. It should be noted that the quantity \G\(vjv Q '^ (see table 1) is temperature 
dependent and that an iteration process is required in the determination of the rotational- 
temperature value. 

The test data are given in figures 8 and 9 as plots of the percent difference between 
the weighted- mean value of the rotational temperatures and the gas temperature as a func- 
tion of the gas temperature. The gas temperature is the wall temperature of the chamber 
corrected in accordance with figure 3. Two points should be noted. The first point is 
that the general trend of the data indicates a uniform variation in the percent difference 
between the measured weighted-mean value and the gas temperatures with increasing gas 
temperature up to approximately 1000° K. There are two data points in figure 8 at 
approximately 1050° K which do not follow this trend. This is attributed to the fact that 
equation (1) as originally derived is assumed to hold for temperatures up to 800° K (ref. 1) 
and, therefore, would not apply rigorously for higher gas temperatures. The second 
point is that the measured rotational temperature is high at ambient temperature and that 
the difference increased with an increase in gas number density. This result indicated 
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Figure 8.- Results of T r tests for temperatures of 300° to 1000° K. Gas number density of 5 x 10 15 cm" 3 . 
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Figure 9.- Results of T r tests for temperatures of 300° to 900° K. Gas number density of 1 x io 16 cm' 3 . 

that the measured rotational temperature was dependent on gas number density. Several 
investigators (refs. 1 and 4) have pointed out that T r measurements made in a static gas 
were higher than the gas temperature. 

Because of the preceding observations and results, tests were conducted to deter- 
mine measured T r dependence on the gas number density. All tests were conducted at 
ambient temperature which ranged from 287° to 303° K with an average of 291° K. The 
measured T r was compared with the corresponding ambient temperature and plotted 
against the normalized number density as shown in figure 10. An analysis of the results 
provided the following empirical equation: 

JL5.7 

AT = log io(i^) 

where AT is the difference between the gas temperature and the measured rotational 
temperature at the point of observation dependent on the number density n, and n Q is 
the reference gas number density arbitrarily selected as 3 X 10^4 cm~3 because the 
results appear to intersect at this value for a unity ratio of T r /T. The number density 
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Figure 10.- Dependence of T r on gas number density. n 0 = 3 x l0 i4 cm' 3 . 

n is limited to the range, 3 X 10^4 cm _ 3 s n ^ 3 X 10^6 cm"^. On the basis of experi- 
mental observations, this increase in temperature is assumed to be a function of number 
density only and to be independent of gas temperature. Therefore, as the gas tempera- 
ture increases, the change in T r due to number density would become relatively small. 

Equation (2) is assumed to be applicable to only the physical conditions of the tests 
reported herein. That is, effects of beam current, beam potential, and the distance of 
the observation point from the exit aperture of the electron gun were not included. All 
observations were made at a point 30 cm from the gun exit aperture, a beam current of 
1100 juA, and a beam potential of 28.5 kV. The greater the distance from the aperture the 
greater will be the beam spreading for a given gas number density. If the gas heating 
should be dependent on the beam electron current density, spreading will be an important 
factor. It should be pointed out that changing the beam current by a factor of 2 did not 
change the measured T r value within experimental precision. Also, the temperature of 
the gas at the point of observation is a function of a temperature gradient and, therefore, 
beam spreading again may be a factor because it could alter the gradient. However, 
until the exact molecular energy exchange mechanism or process is determined, tests 
should be conducted for each test configuration. The molecular heating mechanism must 
be determined in order to know conclusively the effects on the rotational temperature 
measurements in a flowing gas. 

Tests similar to those of reference 4 were conducted with a dove prism. The dove 
prism was used to rotate the beam image so that it was perpendicular to the spectrometer 
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entrance slit. All other tests were conducted with the beam parallel to the entrance slit. 
The purpose of these tests was to obtain some quantitative information about the effects 
of secondary electrons by imaging the beam perpendicular to the entrance slit. Thereby, 
the fluorescence due to primary electrons was reduced and the fluorescence from the halo 
region about the beam center which is attributed to secondary electrons was increased. 
Within the experimental precision of this work, no difference was noted between T r 
measurements taken with the beam image perpendicular and those taken with the beam 
image parallel to the entrance slit. This result did not permit any definite conclusions to 
be drawn about the effects of secondary electrons. 

Results and Discussion 

Corrections in accordance with equation (2) were applied to the measured rotational 
temperatures. The corrected results are shown in figure 11 which is a composite plot of 
data shown in figures 8 and 9. The resultant scatter of data has been reduced from +10 
to -4 percent (figs. 8 and 9) to +2 to -6 percent (fig. 11). Also, the data are approximately 
grouped about a common curve. A slight difference in the data still exists for the two 
different densities. This difference can be attributed to the inaccuracy of equation (2) 
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Figure 11.- Results of T r tests corrected for gas- number-density effect for temperatures of 300° to 1000° K. 



because of the experimental uncertainty associated with the data on which its derivation 
is based. The corrected results still indicate a slightly increasing difference between 
gas temperature and the measured rotational temperature. An analysis of equation (1) 
does not reveal any apparent sources of this difference. There are two possible causes. 
The first is possible experimental error. Careful analyses of all experimental apparatus 
have been made, and reasonable confidence has been established. The second possible 
source is found in the formulation of the theory leading to equation (1). The excitation- 
transition process is complex and no direct account of the effects of secondary -electron 
excitation processes has been made. Also, it was assumed that the rotational levels of 
the ground molecular state are in thermal equilibrium. Because a tempera- 

ture difference does exist between the gas temperature and the measured rotational tem- 
perature, there will exist a population distribution of the rotational levels at the point of 
observation which is at least slightly non-Boltzmann (ref. 6). A slight deviation from a 
Boltzmann distribution would not be detected within the obtainable experimental precision 
but would affect the measured rotational-temperature value. 

Additional T r measurements were made in low-density nitrogen gas. The results 
of these measurements are shown in figure 11 for comparison with those measurements 
made from the nitrogen constituent of low- density air. There does not appear to be any 
significant difference between those measurements performed in air and those performed 
in nitrogen. This result is what would be expected because nitrogen is the primary con- 
stituent of air and collision quenching of the fluorescence of the first negative system of 
nitrogen by other constituents of air would not be a problem unless it was selective and 
greatly disturbed the relative molecular number density population of the excited state of 
nitrogen. 

Rotational-temperature measurements were also made from the data obtained from 
the (0,1) band of N 2 - The data reduction was performed in the identical manner described 
in the preceding sections. The values for used in equation (1) are given in 

table 2. The results are plotted in figure 11 for comparison with other results. No 
apparent difference exists between T r measurements made from the (0,0) band data and 
those made from the (0,1) band data. 

A relative precision of ±5 percent was estimated for these measurements. The 
factors affecting the relative precision of the measurements were the uncertainties in gas 
number density, gas temperature, and rotational temperature values obtained from the 
graphic solutions. A value of ±2 percent was calculated for the effects of density and gas 
temperature, and ±3 percent was estimated for the uncertainty of the measured tempera- 
ture based on standard error values obtained from a large number of independent 
measurements. 


VIBRATIONAL- TEMPERATURE MEASUREMENTS 


Tests 

The (0,1) and (1,2) bands of N 2 were selected for the vibrational-temperature 
measurements because there is no overlapping of second positive systems of N 2 at the 
test temperatures and because these bands are sufficiently intense to permit reasonably 
accurate intensity measurements. The test chamber was set to the desired temperature 
and a spectral trace, similar to that shown in figure 12, was obtained with a relatively wide 
spectrometer entrance slit. A slit width of 250 microns was used for these tests because 
it permitted sufficient instrument spectral resolution to distinguish the N 2 bands and 
provided a strong light level at the spectrometer photomultiplier cathode. The relative 
intensity of each band was determined by mechanically integrating the total area enclosed 
by each band envelope. For rotational temperatures greater than 400° K, it was neces- 
sary to correct the measured areas of each band because of the overlapping of the (0,1) 
band onto the (1,2) band. The relative intensity value, as determined from the spectral 
traces, was applied to the curve shown in figure 13, which is a plot of intensity ratio of 
the (0,1) band to the (1,2) band as a function of the vibrational temperature. This curve 
was calculated from the following equation: 
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The derivation of this equation and a general discussion of the theory leading to its deri- 
vation are given in the appendix. Tests were conducted for gas temperatures of approxi- 
mately 300° K, 700° K, and 1100° K. The measured intensity ratios of the (0,1) band to 
(1,2) band are plotted against gas vibrational temperature in figure 13. The intensity 
ratios were corrected for variation in the spectral response of the detector which was 
determined from a calibration made with a standard lamp. The correction factor was 
0.97. 


Results and Discussion 

Normally, when the electron beam technique is used, the vibrational temperature is 
not known. Therefore, the measured intensity ratio is plotted as a point on the curve in 
figure 13, and the corresponding temperature is determined. The precision of the 
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Figure 12.- Typical unresolved spectrometer trace of (0,1) and (1,2) bands of nJ. T ~ 300° K. 
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Figure 13.- Vibrational temperature as a function of intensity ratio of (0,1) band to (1,2) band of Nj. 
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measured temperature is determined by the precision of the measured intensity ratio. 

The main factor affecting this precision is the ability to account properly for the over 
lap of the (0,1) band onto the (1,2) band. For rotational temperatures below 400° K, 
accounting for the overlap is not a significant problem, but the limitation is the precision 
of the mechanics of measuring the area under the spectral envelope of each band. It is 
interesting to note that in hypersonic test facilities in which this technique can be applied, 
a nonequilibrium condition exists between rotational and vibrational energy states. This 
nonequilibrium condition is such that the rotational temperature is less than 400° K, but 
the vibrational temperatures are much greater. Therefore, a precise determination of 
intensity ratios can be made. 

Another factor enters into the uncertainty of a vibrational temperature determined 
by this method. This factor is the accuracy of the transition probabilities used in deter- 
mining the curve in figure 13 from equation (3). The accuracy of the transition probabili- 
ties is usually taken to be between 10 and 20 percent; a 15-percent value is taken for this 
work. 

The total uncertainty of vibrational temperatures determined with this technique 
may be as large as 20 percent. The results of these tests show that this uncertainty is 
a fairly good estimate because the largest error was approximately 16 percent. 

DISCUSSION OF RESULTS 

Two significant results were noted in the work reported here. Neither result was 
predicted by the theoretical model developed by Muntz (ref. 1) for measuring rotational 
temperatures. 

The first result is that the measured rotational temperature is a function of the gas 
number density. A theoretical explanation for this result is not known but it is apparent 
that there is an energy increase, or redistribution, in the rotational energy states of the 
molecules as the gas number density is increased. This increase or redistribution of 
energy could be a result of the excitation process. It is questionable whether this effect 
is due to the fast primary electrons, the lower energy secondaries, or both. However, 
Culp and Stair (ref. 7) recently performed experiments with low-energy (19- to 300-eV) 
electrons, which indicate that low-energy secondaries could be responsible. A portion of 
the secondary electrons is believed to lie in the energy range, =20 eV (ref. 8). The 
experimental results of Culp and Stair showed that the measured rotational temperature 
was a nonmonotonic function of the electron energy in the energy range of their 
investigation. 

The second result is that the difference between the gas temperature and rotational 
temperature increases with an increase in gas temperature. The general trend of this 
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difference is indicated in figure 11 which shows that the rotational temperatures are lower 
than the gas temperatures. In order to obtain this result, it is necessary for the popula- 
tion of the lower rotational quantum states of N 2 X 2 2 g to be enhanced or for the popula- 
tion of the upper quantum states to be decreased by some mechanism. The results of 
reference 4 appear to indicate a continuation of this trend in the low-temperature region, 
^300° K, in which the rotational temperature is increasingly greater than the gas tempera- 
ture for decreasing gas temperatures. This result is important because it indicates that 
some significant contribution to the excitation process is not accounted for in the present 
theory. It remains to be determined whether this result is due to a second-order effect 
which could be attributed to excitation by the fast primary electrons, to excitation by the 
secondary electrons, or to some other process. 

Despite the differences noted in the preceding discussion, the present theory can 
be used for rotational-temperature measurements over the temperature range investi- 
gated with the expectation that the results will be within ±10 percent of the true value. It 
is also believed that these results are pertinent to measurements performed in a flowing 
gas. To establish the relation between measurements in a static gas and those in a 
flowing gas, additional theoretical and experimental investigations are required. There- 
fore, it is essential to establish a theoretical description of the observed results and, on 
this basis, to ascertain the effects in a flowing gas. 

CONCLUDING REMARKS 

Laboratory measurements of rotational and vibrational temperatures have been 
performed over a range from 300° to 1100° K in static low-density air with an electron 
beam probe. The tests were conducted in a unique test chamber which permitted the 
selection of various combinations of test pressures and temperatures. 

Rotational-temperature measurements were made from the (0,0) band of the first 
negative system of the nitrogen constituent of air. Some tests were also performed using 
the (0,1) band of the first negative system of nitrogen. No substantial differences were 
noted between the results of the (0,1) band tests and those of the (0,0) band tests. 

The accuracy of the measured rotational temperatures was found to vary with gas 
temperature and gas number density. A difference existed between the measured rota- 
tional temperature and the gas temperature. This difference increased with increasing 
gas temperature. With the inclusion of an experimentally determined correction for 
the gas-number-density effect, the rotational-temperature measurements were made 
with an accuracy that varied from approximately +2 to -6 percent with increasing gas 
temperature. 
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Vibrational-temperature measurements were made by using the intensity ratios of 
the (0,1) and (1,2) bands of the first negative system of nitrogen. The measurements were 
found to be within ±20 percent of the gas temperature. 

Langley Research Center, 

National Aeronautics and Space Administration, 

Langley Station, Hampton, Va., November 1, 1967, 

125-24-03-22-23. 
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THEORY 

General 

The test gas used in this investigation was air, and the primary sources of visible 
and near -ultraviolet radiation were the first negative and second positive systems of 
nitrogen. The excitation and emission path for the first negative system is illustrated by 
the energy-level diagram in figure 14. A high-energy electron, designated as a primary 
electron, is emitted by a source and has an inelastic collision with a ground-state nitro- 
gen molecule, ^X^Sg. The molecule is excited to the excited ionized electronic state, 

NoB^£ u , from which it spontaneously radiates and drops into the ground energy state of 
" + 2 + 

the ion, N^X Eg. The intensity and spectral distribution of the spontaneously emitted 

I 

I 

1 



I 

I/" 

v'f = 0 — K| 

N 2 X'2g 

Figure 14.- Partial energy-level diagram of nitrogen. 
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radiation reflects the rotational and vibrational characteristics of the molecules that were 
in the ^X^Sg state. 

Rotational and vibrational temperatures may be obtained from an application of the 
intensity of emission equation 


W - ^n^^nm^nm 


(4) 


where N n is the number density population of the initial level of transition, he i-' rml is 
the energy of the emitted radiation as a result of the transition between states n and 
m, r nm is the wave number of the emitted radiation, and A nm is the transition proba- 
bility of emission. All terms of equation (4) except N n are constants or are dependent 
only on the particular transition involved. Application of equation (4) requires the deter- 
mination of the population and its distribution in the initial electronic state, N^B^Su- 

On the basis of the arguments presented in reference 1, it is assumed that 
is populated primarily by neutral molecules in ^X^Zg. Population contributions from 
other possible origins are neglected. Therefore, the population of the rotational energy 
states of N 2 B 2 2 + is determined by the excitation-transition process of N 2 X-*-Sg 
molecules. 

The excitation-transition process may be described through a Born-Oppenheimer 
approximation of the molecular wave function 


^evJAM “ ^e( r i’ r n) r n ^v( r n)^jAM^’^) 


(5) 


where \p e is the electronic wave function with the ith electronic coordinate r^ refer- 
enced to the molecular axis, i p v is the vibrational wave function with nuclei separation 
r n , and is the rotational wave function which is a function of the Euler angles 

The Euler angles relate the molecular coordinate system to the coordinate sys- 
tem of the fixed point of observation. Also, it is assumed that the interaction of the pri- 
mary electron with the orbital electrons can be described by a coulombic potential. 
Therefore, the following matrix element of this interaction may be used to describe the 
excitation: 


{e D e s B 2 sJv , J , A , M' 


y el 

Z, r n 


epX 1 sJv”J"A"M"] 


( 6 ) 


In this expression, 
primary electrons, 
sent the initial and 


ep and 8p represent the initial and final wave functions of the 
e s represents the secondary electron, X^-Sg and B 2 2u repre- 
final electronic wave functions, v" and v' represent the vibrational 
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states, and J"a"m" and J’A'M’ represent the rotational states. Of course, ^ 

i 

is the coulombic interaction term in which the quantity r-Q is the distance between the 
high-energy primary electron and the ith orbital electron. 


Because the particular transition of interest results in the removal of one orbital 

electron, for simplicity V will be replaced with -i- where the subscript 2 

Lu r^ r 12 

i + 2 + 

denotes the removed electron. The electronic state N^B 2 U is formed from the elec- 
tronic state N2X*2g by the removal of a a u 2s orbital electron (ref. 9). Therefore, the 
electronic state may be expressed as X-^Sg = B^sJct u 2s. Expression (6) is now rewritten 


as 


epe g B^S^v' k' a' M t 


'12 


,|b 2 sJct u 2s]v"k”a’ 


M 


(7) 


Note that K has been used in place of J. This substitution may be made by considering 
the applicable coupling scheme which is Hund's case (b) (ref. 9) and by suppressing the 
spin angular momentum. 

In order to evaluate this matrix element for high-energy primary electrons, a 
plane-wave approximation is made for the primary- electron wave functions, and the inte- 
gration over the primary-electron coordinates is performed. This integration gives 
(ref. 10) 



_ 4tt iq-r 2 

1 - o 


( 8 ) 


where q” is the momentum transfer and rg is the position vector of the interacting 
orbit electron. A series expansion of e lc *’ r 2 gives 


e iq-F 2 


1 + iq . r 2 - 



(9) 


In a first-order approximation, the first two terms of this expansion are retained. But 
the contribution to the matrix element by the first term in equation (9) is zero because of 
the orthogonality of the initial and final states of the molecule. Therefore, equation (7) is 
given by 



( 10 ) 
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In order to evaluate expression (10) further, the vector F 2 is transformed to the 
coordinates of the molecular axis through the dyadic D(0,x,4>) which relates the molec- 
ular coordinate axis to the fixed coordinate system of the point of observation. Therefore, 
expression (10) is rewritten as 


i4ire 2 q 

~W~ 



A’M , |d(0,x,«)| 


k"a”m 




( 11 ) 


where ro is the position vector of the secondary electron with respect to the molecular 

“ 9— 

i47ie*q 

axis. The square of the absolute value of the term — — is contained in an excitation 

q2 

function Ce and will be suppressed in the following equations. 


Evaluation of expression (11) shows that for A' = A" = 0, AK = ±1 transitions are 
allowed, whereas the AK = 0 transition is forbidden, as assumed in reference 1. Also, 
further analysis of expression (10) in the form of expression (11) by using higher order 
terms of equation (9) indicates that transitions for AK = ±3, AK = ±5, and so forth, are 
also allowed. Of course, the magnitude of the cross sections for the excitation-transitions 
for AK = ±3, AK = ±5, and so forth, will determine the contribution to the population of 
the excited state. This contribution is expected to be small. 


The square of the second matrix element in expression (11) is defined as the band 
strength or the vibrational transition probability and is designated P v ' v "* This band 
strength may be approximated by assuming a mean value of the internuclear separation 
and may then be written as 


Pv'v” = 


e s B2 S + r| 


| jp 2 ^u°u2®]) 2 ((v 1 |v")| 


( 12 ) 


The overlap integral of the vibrational wave functions squared |(v'|v")|^ is the Franck- 
Condon factor, q v » v *». Generally, equation (12) is then expressed as 


2 

Pv’v” = | R ij I Qy'v" (12) 

In order to take into account the variation of internuclear separation, a method of r cen- 
troids (ref. 11) is used where r is the expectation value of the internuclear separation 
r n as determined by the vibrational wave functions. Now, the band strength is given by 

Pv’v” = | R ij( r n)| 2 q v > v " ( 14 ) 
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The rotational -line strength is given by the square of the first matrix element of expres- 
sion (11) summed over M" and M' 


I 

m’m" 



(15) 


This term is the Honl- London factor, which is tabulated in reference 9. For the excita- 
tion transition of interest, which is indicated by the superscript a, the relative rotational- 
line strength may be obtained from the ratio of the line strength to the sum of the line 
strengths 



K 

S K " 


I 

K' 


S K ' 

V 


(16) 


The relative line strength for emission, designated P^, is the same as that for excitation 
except that the summation is over K". 

The preceding work of this section has been based on a plane-wave approximation 
for the high-energy incident primary electron. It is also necessary to consider low- 
energy secondary electrons because they may contribute significantly to the number of 
ionized nitrogen molecules. Several observations may be made without making a detailed 
analysis of the excitation-transition process. For electronic states of the homonuclear 
nitrogen molecule, x!s| and B 2 E* rotational levels have symmetric and antisym- 
metric states under nuclear exchange. The symmetry properties of X 1 !^ and B 2 S U 
are the opposite with respect to even and odd rotational quantum numbers, and only tran- 
sitions between rotational energy levels with the same symmetry are allowed. Therefore, 
the selection rules with respect to AK which apply to high-energy electrons also apply 
to low-energy electrons. The only difficulty arises when the first-order approximation 
is not sufficiently accurate and AK = ±3 or higher AK transitions are required to 
describe secondary-electron — induced transitions. 


Vibrational Temperature 

The intensity of spontaneously emitted radiation as a result of transitions between 
two vibrational energy states is given by 


ly’ v'2 = N v’ hc Vv” A v’v” 


(17) 
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where N v » is the number density population of the initial vibrational level, ^v’v'^ is 
the wave number of the resultant radiation of the transition, and A v » v h is the transition 
probability of spontaneous emission for transitions between the initial vibrational level 
v' and a terminal vibrational level v 2 . 

For this work, it is assumed that N 2 Xls| is the primary source of molecules 
which are excited to as a result of inelastic electron collisions. Therefore, 

N v ’ and the resultant intensity of spontaneous emission are dependent on the. distribution 
and the number density population of the vibrational energy states, Ny'^, of N 2 X 1 S+. 

The steady-state relation between N v t and N v t* as a result of inelastic electron-N 2 
collision is given by 



I 


Nv-’PvV^ 


(18) 


where the product C e 


P, T TT 

r v V-J^ 


describes the excitation and transition process between 


Vl 


and v'. The term C e is an excitation function, Pv’v'j is the vibrational band 
strength described in the preceding section, and R is the depopulation rate. The term 
C e /R will be designated C e henceforth. If it is assumed that the N 2 XlSg vibrational 
energy states are in thermal equilibrium, the number of molecules in a given N v '^ vibra- 


tional state is given by a Boltzmann distribution 

N „ - JJo -Ev’i/kTv 
Nv l Qv 


(19) 


where Q, 


■I 


-E v '^/kT v 


is the "state sum" or partition function, 


E v 'i 


= G o( v l) 


he is 


the characteristic energy of the v'j level, Go( v i) vibrational term, N 0 is the 

steady- state population of N 2 X* 2 +, and T v is the vibrational temperature. Therefore, 
with the substitution of equation (19) into equation (18), the dependence of N v » on the 
T v of N 2 XlEg is established, and the relation between the vibrational temperature of 
N 2 Xl 2 g and the intensity of spontaneously emitted radiation (eq. (16)) can be determined. 


It is not necessary to make absolute-intensity measurements because T v may be 
determined from the ratio of intensities of two vibrational bands. The resultant equation, 
as derived from equations (17), (18), and (19), is 


29 



APPENDIX 


Iy k v 2 _ V 1 

Iv' v" 




*vi A Wi 


i v 2 V - E v , ,7 kT v 

D >e 1/ P v > V "A V > v "i/ v ’ v " 

L v b v l v b v 2 v b v 2 


( 20 ) 


Note that C e , N 0 , and Q v canceled in the preceding ratio because these terms are 
independent of a particular v',v'^ transition. Therefore, with equation (20), the variation 
of the intensity ratio as a function of T v may be calculated. As an example, the calcu- 
lation for the intensity ratio for the (0,1) to (1,2) bands is shown in figure 13. The tran- 
sition probabilities used for this example were obtained from references 12 and 13. 


Rotational Temperature 

The rotational number density, N^’ , of N^B^S* is a function of the excitation- 
transition process, depopulation rate, and number density, Nj^, of ^X-^Sg. If it is 
assumed that rotational states of v'j are in thermal equilibrium, N K '^ is given by 
(ref. 9) 


n K” = 


Ny’j' 

q7 


( 2 K 4 + l)< 


- E Ki/ kT 


(21) 


where 


Qr 


-II 


2K" + 1 



is the rotational partition function, 


E K ’| = F(K'{)hc 


K 


is the characteristic energy of a rotational state, F(K'{j is the rotational term, and T r 
is the rotational temperature. A relation to a temperature has been established through 
the Boltzmann factor which explicitly assumes that thermal equilibrium exists in the 
ground electronic state of the neutral species of N 2 . 


In order to interpret this temperature dependence in the resulting intensity of radia- 
tion, the selection rules for transitions between various rotational energy states are 
applied. As a first approximation, the applicable selection rule for the excitation- 
transition process is taken to be AK = ±1. The AK = ±1 selection rule predicts the 
formation of a P branch and R branch in the rotational fine structure of a vibrational 
band in excitation as well as emission. Figure 15 shows the R and P branches of the 
(0,0) band of nJ. 

The following derivation is similar to that of reference 1 by E. P. Muntz. With the 
formation of the P and R branches, the steady-state population of N^B^S^ is given 
by 
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n K’ = Ce J fc-lF 
v" l 


RR + N Ki+1 P RP 


P ] 


p vV{/ 


( 22 ) 


where Ppp and Ppp are the relative rotational- line strengths of absorption, pre- 
viously described, for the P and R branches. 

With the determination of Nr’ , the intensity of emitted radiation may be calculated 
as a function of the rotational temperature T r . Before the calculation may be accom- 
plished, it is necessary to set up an expression for the emission transition probability 

Ak’k’ 2 as 

AK’Kg = X^PvV'Pr (23) 


p 

where X is a constant, v is the wave number of the transition, and Pp^ is the rela- 
tive rotational transition probability for emission. 


P Branch 
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The intensity of emission for a particular R-branch transition is given by 


where 


iR’K’^ = xc e 


*1 


v; (- 


N o ” E v'j A t v 

<3^ V Pv'v'i(A) 


Qr 


K 


,2K' + 1/ V V 2 


■Ett- 1 i /kT r 

(A) = K’e / + (K f + l)e 


" E K , +l/ iT r 


( 24 ) 


(25) 


This equation is simplified by noting that the product 

XC e P vV 'N 0 /Q v = W 


(26) 


where W is a constant for a particular ^v’ jVgJ transition. Also, the equation may be 
put in conventional form (ref. 9) by noting that 2K 1 = K' + k'jj + 1 for R-branch transi- 
tions. Therefore, 


K’ + K '2 + 1 


-E v "/kT 


WiA 


2(2K' + 1) i-t 


Y p v'v’i e 


1 / 


(A) 


Qr 


(27) 


For T v = 800° K, 99 percent of the total population is in the v'j = 0 level. Equa- 
tion (24) may now be written 


^'k'^o 
K f + Kg + 1 


= Yi^jcQe 


-B 0 K’(K'+l)hc/kT r 


(28) 


where Bg is the rotational constant related to the v{ = 0 vibrational level of 
N 2 xis+ 

- E v " /kT v 
WP v v;e V 

Y = (29) 

2Q r I 0 


and 
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[G> 


2B 0 K’hc/kT r 

iV G 


± (K + 1)e -^0(KV 1 )h c /RT r 


2K* + 1 


(30) 


Note that a reference intensity I 0 has been included to permit the measurements of 
relative intensities. The term [g] involves T r and requires a solution of equa- 
tion (21) through a process of iteration. 

Analysis of equation (27) also reveals that only a small error («±1 percent) will be 
introduced by extending the application of equation (28) to the temperature of 1000° K. 
This analysis was performed by calculating an effective value' for rotational constant to 
replace Bq in equation (28) . 

For ease of application, equation (28) is put in the following form: 


^^■K’(K' + 1) + Z = -2.3 log in 
kT r 


iK’K’^o 


(k’ + K 3 + l)j[G](^ 0 )‘ 


(31) 


where v 0 is a reference wave number used to normalize v. From reference 1, the 
v 0 value is chosen for the (3,2) transition. Also, Z is just -2.3 log^gY, which is a 
constant. 
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TABLE 1.- LOG 10 [G](v^ o ) 4 VALUES FOR (0,0) BAND. 


x 

5C.0 

75. C 

100.0 

125.0 

150.0 

1 75.0 

200.0 

225.0 

250.0 

275.0 

390. 0 

325.0 

350.0 

4 CO. 0 

0 

-0.C5C7 

-0.C341 

-0.0258 

-0.0208 

-0.0175 

-0.0151 

-9.0133 

-0.0129 

-0.C1C9 

-0.0099 

-0.0092 

-0.0086 

-0.0080 

-0-0071 

1 

-0.0445 

-0.C312 

-0.0241 

-0.0196 

-0.0166 

-9.0144 

-3.0127 

-0.0114 

-0.0103 

-0.0095 

-0.C088 

-0.0081 

-0.0076 

-9.0068 

2 

-0.0325 

—0. 0257 

-C.02C8 

-0.0174 

-0.0150 

-0.0131 

-9.0117 

-C.C1C5 

-0.CC96 

-0.0080 

-0.0081 

-0.0076 

-0.0071 

-9.0063 

3 

-0.0152 

-0.C178 

-0.0162 

-0.0144 

-0.0127 

-0.C1 14 

-9. 0103 

-C.0C94 

-C.0C86 

-0.0079 

-0.0973 

-0.0968 

-0.0064 

-0.0057 

4 

0.0073 

-0.0C73 

-0.C102 

-0.0104 

-C.C099 

-0.0092 

-0. 0085 

-C. 0079 

-0.C073 

— C. 0068 

-C.0064 

-O.OC60 

-0.0056 

-0.0050 

5 

C. 034 3 

0.CC54 

-0.0027 

-0.C054 

-0.0063 

-0.0065 

-3.0064 

-C.CC61 

-C.CC58 

-C.0C55 

-0.0952 

-0.0049 

-0.0046 

-0.0042 

6 

C. 0653 

C.C203 

0.0060 

C.0004 

-0.0022 

-0.0033 

-9.0038 

-C.C040 

-C.CC40 

-0.0040 

-0.C039 

-0.C037 

-C .0036 

-0.0032 

7 

0.0998 

C.C373 

0.0161 

0.0070 

0.0026 

0.0003 

-9.0010 

-C.CC17 

-C.CC21 

-0.0023 

-0.C024 

-0. 0024 

-0.0923 

-0.0022 

8 

0.1372 

C. C 562 

0.0274 

0.0145 

0.0080 

0.0944 

3.0023 

O.OCIO 

0.0CC2 

-0.0003 

-0.0907 

-0.0009 

-0.0010 

-0.0011 

9 

C. 1 772 

C.C769 

0.0398 

0.0229 

0.0140 

0.0090 

0.0059 

0.0C40 

0.C027 

C.0C18 

0.0012 

0.0008 

0.0005 

0.0002 

10 

0.2192 

0.0991 

0.0533 

0.0320 

0.0206 

0.0140 

0.0099 

0.C072 

C.0C54 

0.0041 

0.0932 

0.0026 

0.0021 

0.C015 

11 

C. 2630 

C. 1 228 

0.0679 

0.0418 

0.0278 

0.0194 

9.0142 

0.0107 

C-0C84 

0.0067 

0.0055 

0.C046 

0.0039 

0.0030 

12 

0.3081 

0.1477 

0.0935 

0.0524 

0. 0354 

0.0253 

9.0188 

0.0145 

0.C1 15 

0.C094 

0.0078 

0.0067 

0.0058 

9.0045 

13 

0.3544 

0.1739 

0. 1000 

0.0638 

0.0437 

0.0316 

0.0239 

0.C186 

C.015C 

0.0124 

0.0104 

0.0089 

0-0078 

9.0062 

14 

0.4016 

C. 2010 

0.1173 

0.0757 

0.0524 

0.0383 

0.0292 

0.0230 

C.C186 

0.0155 

0.0131 

0. Cl 13 

0.0100 

0. 0080 

15 

0.4496 

0.2291 

0.1355 

0.0883 

0.0617 

0.0454 

0.0348 

0.0276 

C.C225 

0.0188 

0.0160 

0.0139 

0.0122 

0.0099 

16 

0-4981 

0.2580 

0.1543 

0.1015 

0.0714 

0.0529 

0.0408 

0.0325 

C.C266 

0.0223 

0.0191 

0.0166 

0.0147 

3.0119 

17 

0.5471 

0.2876 

0.1738 

0.1152 

C.0816 

0.0607 

0.0471 

0.0377 

0.C309 

0.0260 

0.0223 

0.0194 

0-0172 

0.0140 

18 

0.5965 

C.3177 

0.1940 

C.1295 
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